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Following a recent proposal [C. Muschik et. al., Phys. Rev. A 83, 052312 (2011)], engineered 
dissipative processes have been used for the generation of stable entanglement between two macro- 
scopic atomic ensembles at room temperature [H. Krauter et. al., Phys. Rev. Lett. 107, 080503 
(2011)]. This experiment included the preparation of entangled states which are continuously avail- 
able during a time interval of one hour. Here, we present additional material, further-reaching data 
and an extension of the theory developed in [C. Muschik et. al., Phys. Rev. A 83, 052312 (2011)]. 
In particular, we show how the combination of the entangling dissipative mechanism with measure- 
ments can give rise to a substantial improvement of the generated entanglement in the presence of 
noise. 
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I. INTRODUCTION 

In a recent experiment [L ^, a new technique for 
generating extremely robust and long-lived entanglement 
has been demonstrated following a proposal put forward 
in [3 (see also [4 ). By means of reservoir engineering 
[5H26]. entanglement has been produced purely dissi- 
patively. Moreover, it has been shown how engineered 
dissipative processes in combination with continuous 
measurements can be used to create entanglement in a 
steady state. Using this method, entanglement between 
two macroscopic atomic ensembles [27] has been main- 
tained and verified for up to one hour. This extends the 
time intervals during which event-ready entanglement 
of material objects can be provided by several orders 
of magnitude. In this article, we give an extended 
description of this experiment and present additional 
supporting data. A detailed and rigorous derivation 
of theoretical framework of the employed dissipative 
entangling mechanism can be found in [3] . Here, we also 
discuss an extension of this method which has been used 
in the experiment and includes measurements. We show 
how they can be used to improve the purely dissipative 
protocol and explain the basic working principle in detail. 

The coupling of a quantum system to its environment, 
commonly referred to as dissipation, is traditionally con- 
sidered to be a main problem impairing experiments in- 
volving quantum superposition states and the develop- 
ment of quantum technologies. Harnessing dissipative 
processes rather than aiming for eliminating their influ- 
ence is a radically new concept and represents a paradigm 
shift in Quantum Information Science. We show that 
even limited control of the coupling between system and 
environment can enable one to turn a major problem 
into an asset. This change in perspective is not only of 
conceptual interest, but yields also significant practical 
advantages. The protocol discussed here relies on engi- 
neered dissipation. More specifically, the coupling of a 
system with a reservoir is tailored such that the desired 
state is obtained as the steady state of the dissipative 



evolution. This way, the target state is reached indepen- 
dently of the initial conditions. Accordingly, this type 
of protocol does not require the precise initialization of 
the system in a well defined state. The resulting quan- 
tum state can be maintained for long times since it is 
stabilized by the dissipative dynamics. This mechanism 
continuously drives the system into the desired entan- 
gled state, even in the presence of noise sources, which 
limit the coherence time of the quantum system. Thus, 
the use of engineered dissipation enables the realization 
of unlimited entanglement lifetimes, which is not achiev- 
able by traditional methods. 

We consider here two macroscopic atomic ensembles 
at room temperature interacting with freely propagat- 
ing coherent light. Quantum information can be en- 
coded in collective atomic spin states which are unaf- 
fected by the thermal motion of the atoms. This sys- 
tem has been shown to provide an excellent platform for 
quantum memory schemes and the realization of light- 
matter interfaces [28H35] . In the protocol discussed here, 
two atomic ensembles are entangled by virtue of a dis- 
sipative mechanism which is induced by the application 
of a strong driving field. The generated entanglement 
can be accessed at any moment during an extended pe- 
riod of time, which makes it particularly useful for pro- 
tocols, where it is not known in advance when the en- 
tangled state is needed (for example if probabilistic sub- 
routines are involved). If the entangled state is to be 
used, the driving field inducing the entangling mecha- 
nism is switched off before the actual protocol is run. 
Since entanglement is created between two ensembles, 
the resulting atomic state can either be used directly 
or be read out on demand using light-matter interface 
schemes [36l [37] . Another very interesting application is 
the use in continuous protocols, for example in dissipa- 
tive quantum repeater schemes [38 . This type of scheme 
requires continuous entanglement for establishing high- 
quality steady state entanglement over large distances. 

The remainder of the article is organized as follows. In 
Sec. |ll| we summarize the main results and explain the 



key features of the scheme. Sec. HI is concerned with the 




FIG. 1: Setup for dissipative entanglement generation be- 
tween two spatially separated atomic ensembles. A homoge- 
neous magnetic field, which is applied along the x-direction 
defines the quantization axis and leads to a Zeeman split- 
ting of the atomic ground states 1 1) and | ^) . A strong y- 
polarized laser beam (shown in green) couples these states 
off-resonantly to the excited states |e^) and |e^) and to the 
light field in x-direction (shown as wavy lines). The Zeeman 
shift Q of the ground states leads to the emission of photons 
into two sideband modes centered around cjl ±^, where cjl is 
the frequency of the applied laser field. The vacuum modes in 
x-direction provide a common environment for the two atomic 
systems. Due to collective effects, the scattering of photons 
in this direction is enhanced. 



creation of purely dissipative entanglement. We offer an 
intuitive explanation and data supporting this interpre- 
tation. We also compare the basic working mechanism 
to standard approaches and highlight the distinguishing 
features to the dissipative scheme discussed here. There- 
after, a hybrid method is described, where the dissipative 
mechanism is combined with continuous measurements 
on the light field. In Sec. |IV[ we explain how monitoring 
of the scattered photons can lead to an improvement of 
the produced entanglement in the presence of noise. In 

V} we present additional experimental material and 

Vl| concludes the paper. 



II. OVERVIEW AND CENTRAL RESULTS 

In this section, we provide a brief overview to the 
method of dissipative entanglement generation put 
forward in pj [2, and [3^. We explain the main idea, 
introduce the experimental setup and summarize our 
main results. 

Key idea 

As outlined in the introduction, the entangling mecha- 
nism employed here is invoked by reservoir engineering 
and drives the system into a unique inseparable steady 
state. In the absence of other decoherence mechanisms, 
the steady state of the system (corresponding to the re- 



duced density matrix after tracing out the environment) 
can be a pure state. In our case, the reservoir consists 
of a continuum of electromagnetic vacuum modes, 
which provide a common environment for the two 
ensembles. The atomic system can be coupled to this 
environment in a controlled fashion by applying suitable 
laser fields. More specifically, we consider the setup 
shown in Fig. [l] Each atomic ensemble consists of a 
large number N of hydrogen-like atoms with an internal 
level structure with two ground states |t) and \\). Both 
ensembles are driven by a far off-resonant ^-polarized 
coherent field. This strong classical driving field induces 
effective ground state transitions | t) -^ |e|) -^ \ \) 
and II) -^ \e^) -^ |t), which involve the emission of 
f -polarized photons (compare Fig. fl]). This way, the 
classical driving field couples the atomic system to 
the bath of electromagnetic modes in x-polarization. 
The basic entangling mechanism can be understood 
by considering the x-polarized vacuum modes in the 
direction of the laser field with wave- vector Icl and the 
rest of the modes separately. The latter give rise to 
the standard spontaneous emission and represent noise 
processes. The former are shared by both ensembles and 
provide therefore the desired common environment. In 
the setting considered here, emission into the forward 
direction is collectively enhanced for a large optical 
depth d [33 . Hence, these modes can successfully 
compete with all the others and the entangling processes 
happen on a faster time scale than the undesired ones. 

In our case, the target state is a two mode squeezed 
state which is entangled in the collective spin states 
of the two atomic ensembles |^epr)- This state is 
reminiscent of the entangled quantum state introduced 
by Einstein, Podolski and Rosen (EPR) [39 . |^epr) 
is the simultaneous eigenstate with eigenvalue zero 
of two nonlocal operators A and B, ^I^epr) = 0, 
^I^epr) = 0, where 



A 



II Ji — vJii , B = jiJ^ — vJ^ . 



(1) 



J^j are collective spin operators with J =Xli lt)i(i|- 
ji = cosh(r) and v = sinh(r), where r is the so- 
called squeezing parameter [40]. This type of entan- 
gled state is the main working horse for applications 
in quantum information science with atomic ensembles 
and continuous variable systems in general [41]. En- 
tanglement can be verified and quantified by the pa- 
rameter ^ = Sj/(2|(J^)|) = (/i - z^)^ with Sj = 
var(J^j — Jy^u) + var(Jzj — Jz,u)- ^ < 1 certifies the 
creation of an inseparable state [42] [43] . 
The existence of a steady state |^epr) with A|^epr) = 
^I^epr) = can be understood as an interference pro- 
cess where a (x-polarized) photon, which is emitted in 
forward direction with a given frequency could have been 
originated from either of the two ensembles. For a par- 
ticular atomic state, |^epr), these two process interfere 
destructively, such that no photon is scattered into this 
direction. As a consequence, the interaction with the 



laser is switched off and the state remains the same. This 
interpretation can be quantitatively understood in terms 
of a master equation. 

In the scheme reported on here, the system-reservoir cou- 
pling is engineered such that it gives rise to a dissipative 
dynamics which is governed by 



a) 



Ensemble I Ensemble II 

B^ie'd 



dt 



P = ^ent{p) (2) 

(X (ApA'^-A'^Ap+H.C.) + (BpB'^-B'^Bp+H.C.) . 



It can be shown ^ that Pepr = |^epr)(^epr| is the 
unique steady state of this evolution. In order to obtain 
a two mode squeezed state, the system is coupled to 
two reservoirs, which give rise to the jump operators 
A and B respectively. To this end, the ensembles are 
placed in a homogeneous magnetic field which causes 
a Zeeman splitting of the atomic ground states ft. 
Due to the different ground state energies, photons are 
scattered into two different frequency bands centered 
around cjl ± ^, which we will refer to as upper (blue) 
and lower (red) sideband respectively. For a sufficiently 
large separation in frequency space {Q ^ rAtomic^ 
where FAtomic is the largest effective atomic transition 
rate for ground state transitions | t) ^ | |) [3 ), the 
continua of modes in the lower and upper sideband can 
be treated as independent reservoirs (compare IaI). The 
first (second) term in Eq. ^ is due to the interaction 
with the photons in the lower (upper) sideband. Note 
that the jump operators defined in Eq. ^ are nonlocal 
(i.e. involve atomic operators referring to both, the 
first and the second ensemble) and can give therefore 
rise to an entangled steady state. As explained above, 
it originates from the interference between processes 
where a photon is emitted in forward direction by the 
first or the second ensemble. Other processes can be 
included in the form of additional terms in the master 

equation ^p = jCent{p) + >^noise(p), where jCnoise{p), 

summarizes undesired processes such as spontaneous 
emission, collisions and fluctuating magnetic fields. A 
key point lies in the fact that the rate of the entangling 
processes (i2ent(p)) scales with the optical thickness (due 
to collective effects which originate from constructive 
interference involving each individual atom within a 
single ensemble), whereas the rate of the detrimental 
processes does not. Thus, for sufficiently optically thick 
samples, the creation of entanglement in a ste ady s tate is 
possible even in the presence of noise (see Eq. (Bl ) inlBJ). 



Setup and experimental results 

The experiment is carried out using ^^^Cs vapor at room 
temperature. The two-level subsystem is encoded in 
the two outermost hyperfine levels of the 6S1/2 ground 
state within the manifold with total spin F = 4 [44 . We 
identify |t)/ = |F = 4, m^ = 4), |;)/ = |F = 4, tuf = 3) 
and It)// = |F = 4, m^ = -3), |;)/ = |F = 4, m^ = -4) 
(where rap is the magnetic quantum number). The 
atoms are confined in cubic glass cells which are sepa- 
rated by a distance of approximately 0.5m and have a 
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FIG. 2: a) Experimental setup. The horizontally polarized 
probe light shown in green passes through two oppositely ori- 
ented atomic ensembles. The orientation is generated via op- 
tical pumping. The optical pumping beams are depicted in 
blue. Behind the two atomic cells, the *S'2 detector signal is 
processed by the lock-in amplifier (LA) to take measurements 
on the atomic quantum spin components Jy^z in the rotating 
frame, b) shows the optical pumping scheme. The relevant 
laser frequencies and polarizations are indicated in a level 
scheme. 



spatial extent of 2.2cm. Each cell contains 10^^ atoms 
and is equipped with a paraffin-based spin-preserving 
coating. The experimental setup is sketched in Fig. 2a. 
The two ensembles are prepared in oppositely oriented 
coherent spin states (CSS). This is achieved by optically 
pumping the atoms of the ensembles into tuf = ±4 in 
the x-direction respectively. The circularly polarized 
pump lasers are depicted in blue and Fig. 2b shows the 
atomic level structure, indicating laser frequencies and 
polarization. The strong probe beam which is initially 
polarized in ^-direction transverses the atoms in the 
i-direction. Behind the cells, the detection system is set 
up. 

In [1 , two types of results have been obtained. 
Firstly, entanglement has been created purely dissipa- 
tively, demonstrating that this type of processes can be 
harnessed for tasks in Quantum Information Science. In 
this series of experiments, entanglement has been main- 
tained for a time span, which is an order of magnitude 
longer than the time intervals for which entanglement 
could be sustained in this system so far [33]. Entangle- 
ment has been obtained in a quasi steady state rather 
than in a true steady state due to the multi-level struc- 
ture of Cesium. As opposed to the two-level model dis- 
cussed above, atoms can leave the two-level system by 
undergoing transitions to other internal states. The en- 
tangling processes are fast compared to these undesired 
transitions. The desired dynamics with respect to the 
two-level subsystem reaches a steady state, but since it 
is superposed by the slow detrimental processes involving 
atom losses, entanglement disappears. 
In a second series of experiments, strong pump and re- 
pump fields have been applied in order to transfer the 




FIG. 3: Entanglement generated by dissipation a) Entangle- 
ment for different detunings SQ of the magnetic fields. The 
different curves correspond to 6Q — (black), 6Q — 20Hz 
(violet) and SQ = 40Hz (blue), b) Dissipatively generated 
steady state entanglement. Predicted time evolution of ^(t) 
in the presence of strong pump and repump fields for d = 55, 
(blue), d = 100 (grey) and d = 150 (black). The other param- 
eters take values close to the ones used in 1 , as explained in 



atoms back, which left the two-level system. These inco- 
herent fields lead to increased noise contributions, which 
prevents the generation of a purely dissipative steady 
state in this particular setting. However, by combin- 
ing the dissipative mechanism with continuous measure- 
ments on the light field, a true steady state has been 
obtained. In the following, we consider the generation of 
entanglement by dissipation in more detail. 



III. PURELY DISSIPATIVE ENTANGLEMENT 

This section is concerned with the purely dissipative 
generation of entanglement. The basic mechanism is 
explained and further substantiated by additional exper- 
imental data. Moreover, the perspectives for obtaining 
steady state entanglement for multi-level systems is 
discussed. 

Entanglement creation by virtue of interference 
As mentioned earlier, the underlying mechanism can be 
understood as an interference effect in the dissipative 
channel. Due to destructive interference, no photon is 
scattered into the forward direction in the steady state. 
Accordingly, no measurement needs to be performed on 
the light field. Additional evidence that the entangle- 
ment is produced by the collective dissipation process is 



provided by its dependence on the dephasing between 
spin coherences of the two ensembles. To demonstrate 
this dependence, we introduce a detuning SQ = Qj — Qjj 
of the Larmor frequencies of the two ensembles by tuning 
the bias magnetic fields. Fig. [3^ shows that already at 
6Q = 20Hz entanglement disappears. This can be under- 
stood as a consequence of the "which way" information 
(or "which ensemble" information in this case) provided 
by the distinguishability of the photons scattered from 
the two ensembles into the sideband modes cjl ± ^i and 
cjl ± ^11 . Hence the atomic samples do not share the 
same reservoir any more and entanglement disappears. 
Note that the effect demonstrated in the experiment 
can not be understood in terms of photons which are 
spontaneously emitted by the first ensemble interact- 
ing with the second one. Due to the large detuning 
(A = 850MHz), the interaction of a single photons 
emitted by the first ensemble with atoms in the second 
one is negligible in the parameter regime accessible 
in the experiment. Here, atomic transitions are to a 
very good approximation only induced by the classical 
driving field (accordingly, the input-output relations for 
atoms in the second sample are identical with those in 
the absence of the first ensemble). 

Comparison of dissipative entanglement generation 
with other methods 

There exists a large variety of methods for creating 
entanglement between two quantum systems. In 
particular, several methods have been devised and 
demonstrated for entangling atomic ensembles. Even 
though some schemes, which have been experimentally 
implemented share similarities with the one described 
here, they are fundamentally different. In the following, 
we explain this in detail and highlight the features of 
the method realized here in comparison with previously 
demonstrated ones. 

In standard approaches, which are based on a coherent 
interaction followed by a measurement ^H |33l |43l [45] - 
[53], two atomic ensembles A, and B are prepared in 
specific pure states |a)A, and |6)b. An additional system, 
E, which typically corresponds to certain modes of the 
electromagnetic field, is also initialized in a specific 
state, for example the vacuum |0)e. For appropriately 
chosen external parameters such as the frequency and 
polarization of applied laser fields, the interaction of 
system E with A and B gives rise to an entangled state 
1^) = U\a,b)A,B\0)E- If system E is measured, e.g. 
using a beam splitter and single-photon detectors, or by 
means of homodyne detection, the state of systems A 
and B is projected onto an entangled state, |^(e))A,B- 
This state depends on the outcome of the measurement, 
e. If no measurement is performed (which corresponds 
to averaging with respect to the possible measurement 
outcomes [54]) the resulting state is not entangled. For 
instance, the DLCZ protocol [45], yields a separable state 
if the photons emitted by the ensembles are not detected. 



Dissipative methods can be described as follows. |a)A, 
|6)b, and |0)e denote again the initial states of systems 
A, B and E respectively. Due to the interaction of system 
E with A and B, the state |^(t)) = U{t)\a,b)A,B\0)E is 
created, where the dependence of the resulting quantum 
state on the time t is explicitly indicated. Under ideal 
conditions, i.e. if systems A and B do not couple to 
other environments, the interaction of A and B with E 
can be engineered such that the atomic system evolves 
towards an entangled state. In contrast to the schemes 
described above, the implementation of this entangling 
dynamics does not require measurements on system E. 
This type of behavior can occur if system E possesses 
an infinite number of degrees of freedom, such that 
a non-unitary dynamics drives the system towards a 
fixed state. Due to this property, E is typically referred 
to as environment and the corresponding interaction 
with systems A and B is referred to as dissipative 
process. Dissipative phenomena of this kind are best 
described by means of master equations. To this end, 
the environment is traced out and an equation for the 
reduced density operator of systems A and B, p, is 
derived as described in [A] In the presence of other 
environments, the dissipation induced by the coupling 
of A and B to system E can still create entanglement 
with a life time, which exceeds the decoherence times 
due to these extra noise sources significantly, if the 
corresponding (uncontrolled) coupling is sufficiently 
weak. Note further, that typically, noise processes can 
be included in the master equation description as it is 
done in the present work. 

In the experiment discussed here, entanglement in- 
duced by dissipation has been observed. In particular, 
in contrast to approaches which have been previously 
implemented, entanglement is obtained without using 
measurements on the quantum state of the environ- 
ment [55 . Furthermore, systems A and B remain 
entangled for 40ms. This entanglement life- time is at 
least by a factor 16 longer than the decoherence time in- 
duced by other noise sources. It has been experimentally 
verified that in the absence of the dissipative process, 
the measured entanglement life time is limited to 2.5ms 
due to the remaining noise sources such as collisions or 
inhomogeneities of the applied magnetic fields. 
Dissipative methods exhibit another distinctive feature, 
which is present for an ideal two-level system (compare 
[3 ), but not in the multilevel description of the experi- 
ment. For long times t ^ oo, systems A and B decouple 
from the environment E, |^(t)) -^ \^)A,B\E{a,b,t))E 
under ideal conditions, i.e. in the absence of additional 
noise sources. Remarkably, the desired state \^)a,b is 
reached irrespective of the initial state of systems A and 
B which can be highly mixed. Moreover, except for an 
initial waiting time, no special timing is required. This 
behavior is again due to the fact that E possesses an 
infinite number of degrees of freedom, which guarantees 
that revival effects are not present. This way, entropy is 



transferred from the system to the environment, which 
drives A and B into a particular steady state, which 
depends only on the engineered coupling. 
Using dissipative methods, a mixed but still entangled 
steady state can be reached even in the presence of 
additional noise sources, as long as the coupling of A and 
B to other environments is sufficiently weak compared 
to the engineered dissipative processes. This opens up 
the possibility to keep systems A and B entangled for 
arbitrarily long times. 

Perspectives for creating purely dissipative steady 
state entanglement in multi-level systems. 
Below, we investigate the possibility of generating purely 
dissipative steady state entanglement in atoms with 
multi-level ground states. The main reason, why the 
system in the experiment reported on in [l] does not 
display purely dissipatively generated entanglement in a 
steady state is the depopulation of the relevant two-level 
subsystem due to spontaneous emission which transfers 
the atoms into other Zeeman levels. The depopulation 
of the relevant levels can be avoided by applying strong 
laser fields, which transfer atoms back. However, these 
fields introduce additional decoherence processes which 
inhibit the creation of entanglement in a steady state. 
This problem can be circumvented by increasing the 
optical depth of the atomic ensembles such that the 
entangling dissipative process prevails over the noise 
processes and dominates the dynamics. 
We consider this scenario by including additional (j± po- 
larized pump and repump fields, which induce resonant 
transitions with Amp ± 1 in the first/second ensemble 
in the model. Pump fields drive transitions within the 
manifold of atomic states with F = 4 and repump fields 
transfer atoms from states with F = 3 back to F = 4. 
As explained in [C) we estimate the effect of these fields 
using a simplified model, which has been used in [1 to 
fit the experimental data. Fig. ^jp shows the predicted 
time evolution of entanglement for d = 55; 100; 150 in 
the presence of both, pump and repump fields (in [1 , an 
optical depth of d = 55 was used) . The other parameters 
used in this calculation take values close to the ones 
used in [1 . The theory predicts that under the present 
maximal optical depth (i = 55, the steady state atomic 
variance is just slightly above the separability criterion 
(this has also been confirmed experimentally). However, 
the experimental realization of purely dissipative steady 
state entanglement should be feasible along two possible 
routes. Firstly, atoms possessing two-level electronic 
ground states can be used, for example Ytterbium 
(^'^^Yb) [564 iSIl • In this case, the two-level theory 
formulated here can be directly implemented, avoiding 
additional dynamics which leads to the growth of ^ with 
time. Alternatively, a true dissipatively generated steady 
state using multi-level ground states can be achieved 
for higher optical depths, which can be obtained, for 
example, by placing the atoms inside a low finesse 
optical cavity. 
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of the quadratures 

^IjII = JyJ/Ill \J\{JxJ/Il)\i 

Pi/ii = ^Jzj/ii/y\{Jx,i/ii)\- 

Light propagates in i-direction (see Fig. fl]) and inter- 
acts with the atomic ensembles. We consider here a one- 
dimensional model which includes only light scattered in 
forward direction. Processes corresponding to scatter- 
ing of photons into other directions enter in the form of 
noise. Light is characterized in terms of spatially local- 
ized modes [59l |60] 



y{z) = -^Jduj(^a^e'^^^-^->^H.c}j 



FIG. 4: Steady state entanglement assisted by mea- 
surements, a) Illustration of the interaction of atoms 
and light in terms of spatially localized modes b) Squeezed 
atomic variance var {Pc/s) in the steady state versus the ratio 
7extra/7s in the absence of measurements (dashed line) and if 
the ^-quadrature of the scattered light field is measured (full 
line) . 7extra and 7s denote the rates of the desired entangling 
processes and the atomic decay respectively. 



In [T], we devised and implemented an alternative 
approach, which enables the creation of entanglement 
which persists for arbitrarily long times. This alterna- 
tive approach combines the dissipative mechanism with 
continuous measurements as explained below. 



where c is the speed of light, b and cjl are the bandwidth 
and central frequency of the applied laser field. The op- 
erators for spatially localized field modes y{z) and q{z) 
obey the canonical commutation relation [y{z)^q{z')] = 
ic5}){z — z')^ where the deltafunction has a width of the 
order c/h. The spatial argument z refers to the dis- 
tance along the propagation direction z. Atoms and 
light interact according to H = Ha + i^L + ^int ^ where 



iJA = f(X| + P|)-f(X|,. 



Pjj) describes the Zee- 



man splitting of the atomic ground states. Due to the 
applied magnetic field, the transverse components of the 
collective spin described by Xj/n and Pj/n precess at the 
Larmor frequency Q. Hi, represents the free propagation 
of light -^viz) = i[Hi,,y{z)] = —c^y{z). The interaction 
of the light field with two pointlike atomic ensembles [61] 
located at 2; = and z = R is given by 



IV. DISSIPATIVE ENTANGLEMENT ASSISTED 
BY MEASUREMENTS 

In the following, we show by means of a simple model 
how measurements on the light field can improve the 
generation of entanglement under the dissipative dynam- 
ics described above in the presence of noise sources. In 
the basic model employed here to illustrate the relevant 
effects, the Holstein-Primakoff approximation [58] is 
used to describe the atomic system and noise is included 
in the form of decay of the transverse spin components 
Jy and Jz at a rate 7extra (a more detailed discussion 
is to be published elsewhere). In this subsection, the 
dissipative generation of entanglement assisted by 
measurements is explained in terms of input-output 
relations, since this approach is more illustrative than 
the master equation formalism employed in Sec. |ll] Both 
descriptions are equivalent and yield the same results. 

For large, strongly polarized atomic ensembles, col- 
lective spins can be described by bosonic modes in terms 



H]r^t = 



/27, ( -Xjy{0) + ZPjq{0) 




^27, ( -Xiiy{R) + ZPiiq{R) ] , 



where Z = fi -\- u and 1/Z = /j, — v. We assume in 
the following, that the distance R between the ensembles 
can be neglected, which is justified for FAtomic^ <^ c, 
where FAtomic is the effective rate at which transitions 
between ground states occur, and L'^ki, ^ R, where L 
is the spatial extend of an atomic ensemble and kL is 
the wave vector of the driving field. The former is a 
necessary condition to neglect retardation effects, while 
the latter is used in the calculation of averaged emission 
rates for fast moving atoms (compare [A| and [3 ). Both 
conditions are well fulfilled for the parameter regime of 
the experiment, where the effective decay rates are of 
the order of few ms and R is on the order of a meter 
(L is about 2cm and kL on the order of 10^m~^). The 
resulting Heisenberg equations can be solved by changing 
to a coordinate system rotating at the Larmor frequency 



Q and performing the variable transformation ( = ct — z 
such that y(C, t) = y{ct — C^t)- We introduce exponential 
cos(Qt) and sin(r^t) modulated light functions 



2/c± 



^f±Jo 



,±7at 



cos(m)y(C,0), 



9i"± = ^/ e^-'^'cosinmCO), 



(3) 



with A/+ 



_ 2^77 



and A/"- = 



?/^"± and 



g'g"^ are defined analogously. We assume that the Lar- 
mor precession is fast compared to the atomic evolution 
{QT ^ 1), such that the operators describing the sin 
and cos modulated light modes are canonical and inde- 
pendent, [^c,±,^c,±] = [ys,±^qs,±] = i and [^c,±,^s,±] = 
[ys,±,qc,±] =0. This yields 
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II) 



and 



-4- (X/ — Xji) was used. As 



where X^ = ^ (X/ 

a next step, we include continuous measurements on the 
light field and consider the corresponding time evolution 
of the atomic state in the Schrodinger picture. The con- 
tinuous interaction and measurement process shown in 
Fig. [l] is illustrated schematically in Fig. |4^ in a dis- 
cretized way. Spatially localized light modes correspond 
here to infinitesimally short pulses of duration r ^ 1/b 
(where b is the bandwidth of the incident laser field as 
explained above), which interact successively with the 
atomic system. Each of these spatially localized light 
modes is initially in the vacuum state, such that the 
quantum state at time t = nr is given by |^(t))^|0)L,n+i, 
where \'^{t))A denotes the atomic state at time t. Then 
atoms and light are subject to an entangling interaction 
resulting in the quantum state e~*^^|^(t))A|0)L,n+i- Fi- 
nally, the ^/-quadrature of the light field is measured, 
yielding the measurement outcome i/n such that |^(t + 

^))^ = -/Wr^ L,n+l(^n|e"'^10)L,n+l|^W)A, whcrC 

P{yn) is the probability to obtain the result z/n- The 
resulting expression can be expanded up to first order 
in the parameter r yielding a differential equation for 
the time evolution of the atomic system. The atomic 
state obtained after the measurement depends on the 
measurement outcome z/n- We consider here Gaussian 
quantum states, i.e. states with a Gaussian Wigner func- 
tion. These states are fully characterized by their first 
and second moments, - prominent examples include co- 
herent as well as two mode squeezed states. The Gaus- 
sian character of a state is preserved under the evolu- 
tion of Hamiltonians, which are at most quadratic in 



the system operators and under Gaussian measurements 
such as homodyne detection. Since we consider all in- 
teractions and measurements to be of this kind and all 
states to be Gaussian, the entanglement of the resulting 
state is completely determined by the atomic variance 
var(Pc/5), which does not depend on i/n [62. There- 
fore the resulting entanglement is independent of the 
measurement outcome. If the measurement results are 
traced out (p(t+r) = ^^^ Mnp{t)Ml, where M^ =L,n+i 

{yn\^~^^^ \^) L,n-\-i) 1 and the resulting expression is evalu- 
ated to first order in r, the master equation used in Sec.[ll| 
and [A] is recovered, if the spin operators are replaced by 
creation and annihilation operators within the Holstein- 
Primakoff approximation [58 . The whole process can be 
conveniently described by means of the Gaussian formal- 
ism, where atomic states are expressed in terms of their 
covariance matrix Fc, and displacement vector D, which 
display the second moments (variances and covariances) 
and first moments (mean values) of the system respec- 
tively. In particular, this formalism allows one to easily 
calculate the variances of atomic quadratures after the 
Gaussian measurement of the ^/-quadrature of the light 
field, var (^c/s) ^ond ^^ ^^^ ^^^ of each time step depend- 
ing on the variance prior to the measurement 



var (P^ 



c/s, 



cond 
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var(Pe/J ^— ^ — - — ^— ^ , (5) 

' 4var(^c/s) 



where ?/c/s and g^/s refer here to the localized light mode 
interacting with the ensemble in the n*^ time step and 
75r <C 1 is assumed. This way, a differential equation 
for the squeezed atomic variances is derived. In the ideal 



case. 



var(P,/,)^^^^ it+T) = var(Pe/,;^^^^ 



(1 - var (Pe/,) (t)Z2) 7,T, 



which yields 



var 



(^'=/«)cond(*) 



e-7.*(var(P,/,)(0)) \z'^(\- 



whereas in the absence of measurements, 



var (P,/,) (t) = e-^^' var (Pc/.)(0) + ^ (l - e 
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Both time evolutions result in a steady state with 
var {Pc/s) = 1/Z'^ = (/i — vY ^ since atoms and light 
decouple for t ^ oo. Accordingly, the steady state en- 
tanglement can not be improved by means of measure- 
ments on the light field in the ideal case. The situation is 
different in the presence of noise sources, which prevent 
the decoupling of atoms and light. In this case, residual 
atom-light correlations persist in the steady state and 
measurements on the light field can be used to improve 
the entanglement. Here, we illustrate this effect by in- 
cluding atomic transverse decay at a rate 7extra [63 . If 



the ^/-quadrature of the scattered hght field is measured, 
one obtains 



Scond,C 




2Z2 7. 
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Fig. [4]3 shows that for 7extra > 0, the steady state en- 
tanglement described by this equation is higher than the 
steady state if no measurements are performed, which is 
given by 



^o 



^2 7s 



■ Text r a 



7s + Textra 



Note that in principle, all measurement results y{t) ob- 
tained during the continuous measurement procedure 
could be used to perform feedback operations which sta- 
bilize the atomic state at a certain position in phase 
space. However, this is not necessary here, since the 
atomic quantum state at time t depends only on the 
recent history of the measurements, i.e. on y{f) for 
tss ^ t' < t^ where tgs is the time it takes to reach 
the steady state. For the dissipative processes consid- 
ered here, the atomic state p{t) is memoryless regarding 
events which occurred in a time interval longer than tss- 
Thus, only measurement results obtained during a fixed 
time interval tss , which is independent of t are needed to 
localize the atomic state in phase space. 



V. EXPERIMENTAL ATOMIC STATE 
RECONSTRUCTION 

The experimental setup is depicted in Fig. [2] and de- 
scribed in the caption and the surrounding text. This 
section is focussed on the characterization of the atomic 
state via measurements on light which has interacted 
with the atomic ensembles. 

The light observable of interest is the Stokes operator 
5*2 given by the difference of the number of photons po- 
larized in ±45°. In the specific setting discussed here, 
the probe beam is strongly polarized in ^-direction and 
the f -polarization represents the polarization mode of in- 
terest for measurement as discussed in Sec. IIVI Then the 
Stokes operator 5*2 ~ vW^ • V^ where (/) is the photon 
fiux. S2 can be measured with polarization homodyning 
techniques as depicted in Fig. [2] The light beam is send 
through a halfwave plate and a polarizing beam split- 
ter (PBS). The signals from the detectors situated in the 
output ports of the PBS are subtracted and the differ- 
ence signal is analyzed at the Larmor frequency Q with 
a lock-in amplifier, since we are interested in detecting a 
signal from the spins in the rotating frame. Additionally 
the measurement outcome can be weighted with suitable 
mode functions f{t) and after suitable normalization we 
are thus able to measure the light observables yc,s- oi" 
yc,s-\- defined in Eq. (|3|). 



A. Atomic state reconstruction 

Now, the variances of the collective atomic operators 
P'J" = l/V2{Pf + P}"}) and Pf = 1/ V2(Xf - Xf}) can 
be found from measurements on the transmitted light. 
From the variances of the outgoing operators yc,s-^ using 
the input-output relations given in Eq. Q, it follows: 



var(P») 



-^(var(yr)-^i„(l-^)), 



Z2' 



var(Pr) = -^(var(t/r)-^i(l 



Z2 



)), 



(6) 



where af^ is the shot noise of light, assuming that the in- 
coming light is in a coherent state. The coupling constant 
is defined as />: = Z^/l — e~^^^^. The normalized EPR 
variance of atomic noise is ^ = var(P^") + var(Pj") = 
var(Pf + P|5)/2 + var(Xf - X}})/2. 



B. Including decay and detection efficiency 

To include the decay of the atomic spin due to spon- 
taneous emission, magnetic field instabilities, etc. into 
the input-output equations (Eq. [4| , we assume a decay 
with the rate 7extra towards the CSS. Since the interac- 
tion times which are used to perform the read out are 
short, this is an adequate approximation. The atomic 
part of the input-output equations becomes 

pout^pin .g-7T__^ in , -^e\/l-e-'^^^ ■ F I ^ 

^ C,S ^ C,S ^ r72 ^C,S+^'^ V J- ^ -^p,C/S,+ 'i 



(7) 

with the two-cell noise operators Pi+ = 

l/NFj^e-^^^-'^Fi{t)dt and {F^) = \ and 
e^ = 7extra/7- The relevant light mode is exponentially 
growing with the total decay rate 7 = I/T2 = 7s +7extra- 
The coupling const ant is reduced due t o the decay and 
defined as: n = Z ^ (\ - t^)(\ - e''^^^) . 
Also the equations for the light are adjusted accordingly: 



vZ- 



out 

^c^s — 



- e\/l - e2Z(Fp_ - v/1 - kV^'-F;>+), 



- e^/l^Z(F,_ - ^l-K^IZ-^F,^). (8) 



The non orthogonal exponentially growing and falling 
light modes are now mixed due to the decay. 

The reconstruction equation (Eq. [6| must then be ad- 
justed accordingly: 

var(Pi") = 4(var(y?- ) " C^' ' <« " y^P})). 
var(Xr) = 4(var(yr) " t^' ' <,„ " y^P})). (9) 



with the corrected, reduced couphng constant and 
U'^{h?,T2) and V'^{n^,T2) which can be calculated di- 
rectly from Eq. (Is]) [64 . To extract the atomic noise 
from the light noise measurements, it is therefore only 
necessary to know the coupling constant a^, for which a 
measurement procedure is explained below and the decay 
time T2 which can also easily be measured. 

Additionally, the detection efficiency 77 = 0.84(4) which 
arises from light losses and unperfect detection can be 
included by assuming a beam-splitter with transmission 



C. Measurement of the coupling strength 

The most important experimental parameter for the 
reconstruction is the coupling constant n. It would be 
possible to determine h? ^ by performing noise measure- 
ments on known atomic states, e.g. the CSS or the ther- 
mal atomic state. However, imperfect state preparation 
or additional noise sources can spoil such measurements. 
The approach we implement here is therefore based on 
measurements of mean values as opposed to noise mea- 
surement. The modus operandi is to transfer a coherent 
light state with a known displacement to the atoms and 
then read out the atomic state [341 165] . 

First, a pulse is sent through two oppositely oriented 
atomic samples with a displacement in q]^^^ so in Sz- 
Following EqJ4]and assuming that the atoms possess no 
initial displacement, this leaves the atomic sample with 
a mean value in the X-quadratures. 



a) 



(^cT) = ^(C.V)- 



(10) 



To be able to read out those atomic mean values via a 
measurement on yc^s and thus gain information on /^, we 
apply a 7r/2-pulse to the atomic spin rotating X into P. 
This can be done by adding a magnetic field in the x- 
direction, so that the spins rotate a little faster or slower 
in between the pulses and Jzj -^ Jyj and Jyji -^ Jz,ii^ 
leading to Xc^s -^ Pc.s- Then we send a second light pulse 
for the read out. The outcome of the light measurement 
of the second pulse reveals 



(<"/-) = M^r;) 



^'«:V)- 



(11) 



The coupling strength can be calculated: n^ = , ^'/t ^ . 

The displaced coherent light states are produced with 
the help of an electro optical modulator (EOM) [66 . The 
strongly polarized beam is sent through an EOM whose 
optical axis is slightly tilted compared to the input po- 
larization. Then a DC voltage and a small modulation 
at 322 kHz can be used to rotate a small portion of the 
large polarization component in ^-direction into the x- 
polarization mode. The value of the DC voltage deter- 
mines the phase of the modulation in phase space, so the 
position in the S2-S3 plane. The strength and phase of 
the RF modulation determine the size of the displace- 
ment of the cosine and sine modes. 



Pulse sequence: 
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FIG. 5: In a) the pulse sequence is shown. The probe light 
is turned on and off smoothly to avoid noise contributions at 
322kHz. b) shows n^ for a varied number of atoms and a light 
power P — 5mVK, in c) the corresponding decay constants 
7extra and 7s are shown. In d) the scaling of the normalized 
atomic noise over N is shown. 



For the measurement of {q].^s)^ a | -plate is inserted in 
the detection path to switch to a S^ measurement j67] . 

In Fig. [5] measurements of n^ are shown for different 
numbers of atoms N with a fixed probe power of 5mW 
and a probe duration of 1ms. N is varied by changing 
the temperature of the atomic ensemble. It can be mon- 
itored by sending a weak linearly polarized probe beam 
in the direction of the macroscopic orientation. Due to 
the Faraday effect the light polarization is rotated by 
(X Jx and for the CSS Jx ^ 4 ■ N. The orientation 
of the atomic ensembles can be tested via magneto opti- 
cal resonance spectroscopy [68]. Orientations of 0.997(3) 
are reached regularly in this experimental setup. Fig. (Sj) 
displays how the measured decay rate 7 = I/T2 can be 
decomposed in 7^ and 7extra- Clearly 7^ (x P ■ Of- 



D. PN measurement 

When the coupling constant is known, Eq. ([9| can be 
used to reconstruct the collective atomic operator noise 
from measurements of noise on yc,s- of the outgoing light. 
In Fig. [5]i measurements of the atomic noise in units of 
projection noise (PN) are shown for different N. There 
is a small additional noise component, probably arising 
from technical noise, which scales with the number of 
atoms. Only a small range of Op and thus N is shown. 
For higher N additional classical noise of unknown origin 
is measured, disqualifying higher atomic densities as a 
working point for quantum noise limited measurements 
in this setup. 

To find the time evolution of the atomic noise, while 
probe light is present (as the curves shown in Fig. pp 
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at time t 
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FIG. 6: Pulse sequences for long probe durations. On the 
left a pulse sequence with a long probe pulse is depicted, dis- 
playing the evaluation time interval [t; tprobe] used to find the 
atomic noise at time t. On the right it is shown how a mea- 
surement prior to t can be used to conditionally reduce the 
atomic noise. For measurements like this typically a fast ex- 
ponentially growing mode function is used. 



behavior is rooted in the additional decoherence mech- 
anisms. The atomic state acquires a noise component 
piling up with rate 7extra- The measurements closest in 
time to the atomic state which one wishes to squeeze, 
should be weighted most. 

The actual measurement of var(P°y^*)con is done by 
evaluating ?/°y* in the second time slice [t;t-\- tprobe]- The 
same reconstruction mechanism as for the unconditional 
atomic state reconstruction is used, therefore is the sec- 
ond time slice evaluated with an exponentially falling 
mode function with 7 = 7s + 7extra. Again var(?/°^^_) 
is utilized to establish the value of var(P?/^^*) at time t. 



Accordingly, var(y°^*)con 



var(y°^* 



V5--«*^c/f') can be 
used to find var(P°y^*)con with the same reconstruction 
procedure. The conditional variance can be extracted 
from the two-time correlation functions (^^/° ^y^jl_)'- 



Y^T{ylJl_ 



by optimizing a*, this is assuming that there are no cor- 
relations between the light operators at different points 
in time < y{t)y{t') >= S{t - t'). 



yc,s- is evaluated in the time-interval [t^t-\- tprobe]^ where 
tprobe is the evaluation time with /^(tprobe), to find the 
atomic variances at time t. The corresponding pulse se- 
quence is shown on the left of Fig. [6J 



E. Measurement of the conditional atomic noise 
variance 

The theoretical background for generation of the 
steady state of atoms with the noise reduced by com- 
bining dissipation with the measurement is described in 
Sec. |lVl The experimental procedure for the generation 
of a state with reduced conditional noise variance is de- 
scribed in the following. We wish to squeeze the atomic 
variances at time t by measuring the outgoing light op- 
erators yc and ys for the preceeding time. The time se- 
quence is illustrated on the right of Fig. [6] The probe 
pulse is divided in sections, where measurements in the 
first time interval [0, t] are used for the conditional noise 
reduction of the atomic operators at time t, and the con- 
secutive time slice [t^t-\- tprobe] is used for the reconstruc- 
tion of the atomic state as described above. Here, the 
quantity of interest is the conditionally reduced atomic 
variance at time t: var(P°y^*)bon = var(P°y^* — Q^^^/° ^) = 

^2c/s^) (compare 



Eq. |5|. The superscript "probe" refers to the first time 
slice from to t. To achieve an optimal noise reduction, 
a and the temporal mode function of the probe section 
are optimized. The optimized mode function is exponen- 
tially rising with a rate that is typically faster than I/T2. 
This means that only the last bit of the long preceding 
pulse is used for the conditional noise reduction. This 



VI. CONCLUSIONS 

The method discussed in this article allows one to 
deterministically generate entangled states between two 
atomic ensembles over a macroscopic distance. The de- 
sired quantum state is thereby stabilized by the dissipa- 
tive mechanism, which renders the created entanglement 
robust and long-lived. Purely dissipative entanglement 
generation does not require postselection nor condition- 
ing. 

For two- level systems, the dissipative mechanism ex- 



plained in Sec. |in| can be used to produce steady state 
entanglement, i.e. entanglement which is available per- 
manently. This is even possible in the presence of im- 
perfections and noise sources. The use of atoms with 
a multi-level ground state complicates the generation of 
steady state entanglement since atoms can undergo tran- 
sitions to internal levels, which are not coupled to the 
engineered dynamics in the desired way. Atoms with a 
two- level ground state such as Ytterbium (^'^^Yb) appear 
therefore best suited for the realization of the discussed 
scheme. However, steady state entanglement can also be 
created in ensembles of atoms with a multilevel struc- 
ture if strong pumping fields are applied, which trans- 
fer the atoms (incoherently) back to a specific two-level 
sub-system. However, this strategy requires a high op- 
tical depth of the atomic sample, or the use of a low fi- 
nesse cavity. Therefore, a hybrid solution has been imple- 
mented in [1 , which combines the dissipative mechanism 
As explained in Sec. 



IV 



measure- 



wit h measurements, 
ments on the light field yield information on the atomic 
state and can therefore be used to obtain a quantum state 
of higher purity and accordingly produce higher entan- 
gled states. 
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The mechanism considered here is apphcable in any 
system which can be coupled to two photonic sideband 
modes according through a tunable quadratic interac- 
tion, for example for the generation of entanglement be- 
tween two mechanical oscillators. Even though the lim- 
ited coherence times of quantum systems impose typi- 
cally severe restrictions on the life times of quantum su- 
perposition states, dissipative methods for quantum state 
engineering allow one to produce event-ready quantum 
states for applications in Quantum Information Science. 
Exploring and exploiting all advantages of dissipative ap- 
proaches will require both, devising new protocols which 
are capable of generating and processing steady states 
as well as finding realistic and practical ways of imple- 
menting the required coupling of physical systems to a 
bath. 
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Appendix A: Derivation of the master equation 

The interaction of atoms and light illustrated in Fig. [l] 
can be described by the effective ground state Hamilto- 
nian 

H = Ha + i^L + ^int 1 

where excited states have been eliminated using the fact 
that the detuning |A| is large compared to the Doppler 
width fooppier and atomic decay rates Fatomic- Here and 
in the following, Fatomic denotes the largest effective rate 
for atomic ground states including single particle as well 
as collective rates (see below). Ha = ^{Jx,i — Jx,ii) 
accounts for the Zeeman splitting of the atoms in the 
external magnetic field and Hi, = J dk ujk a^ak is the 
free Hamiltonian of the light field, where a\^ is the an- 
nihilation operator of a photon with wave vector k and 
frequency uJk. In a rotating frame, the interaction Hamil- 
tonian is given by 



N N 



^Ak r. 1 i 






-H.C. 



(Al) 



where Ak specifies the the two orthogonal polarizations 
of the light mode with wave vector k. The first and 



second integral cover narrow bandwidths Auis and 
Aujus centered around the lower and upper sideband 
respectively. The atomic operator cri/u^i = |t)////,i(i| 
refers to a particle in ensemble //// at position r^, 
Ak = ki, — k and kL is the wavevector of the applied 
classical field. ^(k)/i and ^(k)i^ denote the effective 
coupling strengths of the passive (beamsplitter-like) part 
of the interaction and the active (squeezing) component 
of the Hamiltonian respectively. The laser field covers a 
very narrow bandwidth around the central frequency cjl 
and is sufficiently off-resonant such that the interaction 
is well within the dispersive regime and absorption 
effects can be neglected. 



Starting from Hamiltonian (Al), a master equation 



of Findblad form can be derived for the reduced atomic 
density matrix p{t). To this end, the Born Markov 
approximation is used, which is well justified for optical 
frequencies. Since the Farmor splitting exceeds well 
atomic decay rates Q » Fatomic, the two sideband 
modes can be treated as independent baths. The effect 
of atomic motion gives rise to noise terms and can be 
included in the master equation in the form of averaged 
coefficients, where the average in time corresponds to an 
average in space. This is legitimate in the fast motion 
limit, where the time scale set by the average velocity 
of the atoms v is fast compared to the time scale of the 
radiative decay Fatomic^ <C 1. In this case, the emission 
of light can be described independently of the evolution 
of atomic positions. 

Using the definitions 



1 ^ 1 ^ 

7x7 E^^^^- ^77X7 E^^^^- 



A/'• 



1 ^ 1 ^ 



the resulting master equation can be cast in the form 

|p(t) = ^drAp{t)A^+^drBp{t)B^ 

N 



i=l 
1 ^ 



+.. , (A2) 

where F is the single particle decay rate and a short 
hand notation was used. Master equations of Findblad 
form f^p{t) = \ [Ap(t)A^ - A^ Ap(t)) ^H.C. with decay 
rate 7 and jump operator A are abbreviated by the 
expression ^p(t) = ^Ap{t)A^ + ... . d denotes the 
resonant optical depth of one atomic ensemble. Note 
that the entangling terms in the first line are enhanced 
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by a factor (i, such that for sufficiently opticahy thick 
samples, the additional noise terms reflecting thermal 
motion are small compared to the desired contributions. 

Next, additional cooling and heating processes, as 
well as dephasing are included. The full master equation 
is given by 

I Pit) = ldTApit)A^ +^-drBp{t)B^ (A3) 

N 



1 ^ 

1 ^ 



+ ... , 

with cr^^j/ii^i = |t)////,i(t| and cru,i/ii,i = \i)i/ii,i{i\- 
The noise terms proportional to F/i^ and Tu'^ in Eq. (A2) 
have been absorbed in the second and third line of 



Eq. (A3) respectively. The last three lines represent sin- 
gle particle processes. Hence they do not feature a collec- 
tive enhancement factor as the entangling terms in the 
ffist line. As shown in [3], Eq. (A3) includes all terms 



that need to be taken into account. Collective dephasing 
terms, as well as collective contributions due to pump 
and repump fields can be neglected. Similarly, the dis- 
tance R between the two ensembles does not play a role 
for kL ^ R/L'^ and kLL :^ 1, where L is the spatial ex- 
tend of an atomic ensemble. For the experimental setting 
under consideration this condition is fulfilled. 



where d{t) = d N2{t)/N, f = Fcooi + Theat + Tdeph and 
P2{t) = 2{J^{t))/ {N2{t)). FoTt^oo and TVs = TV, 



^J2,o 



N 



f + C^FP2,oo 



Next, the time evolution of the longitudinal spin is con- 
sidered. Eq. (A3) yields 



dt 



{Jx{t))2 = — (Fheat + Tcool) {Jxif))'} 



N2{t) 



\^ cool J- heat j , 



such that for constant particle number N2 = N, 

-'-^ J- cool J- heat 



\Jx)2,oo — 
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and therefore 
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f^dTPl^{^-vf 
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cool 



heat 



(Bl) 
(B2) 



in the steady state. In the limit (i ^ 00, this equation 



reduces to ^2,. 



i^-iy) . 



The variation of A^2(^) and P2{t) is slow compared 
to the evolution of T,j^(t). In the limit where the 
entangled quantum state follows the changing particle 
number and atomic polarization adiabatically, ^2(^) is 
given by 



t u\ ^ ^^2(0) -(Ud(t)rP2(t))t 
^^ ^ 2P2W 

^ r^d{t)rP2{t)\fi-uy 
P2{t)(r^d{t)rP2{t) 



(B3) 



(= 



-(f+d(t)rP2(t))^ 



)■ 



Appendix B: Calculation of entanglement 

Below it is shown how the entanglement measured by 
the quantity ^ = T^j/ {2\{Jx)\)^ where Sj = var(J^j — 
Jy,u) + var(Jzj — Jz,u)^ can be determined in the hmit 
A^ ^ 1, assuming that the number of atoms in the two- 
level system N2 depends on time. For clarity, operators 
referring to the two-level model are labelled with sub- 
script "2". The time derivative of the variance T^j^ is 
calculated using Eq. ( |A3[ ). By applying the decorrela- 
tion approximation {Jy/z{t)Jx{t))2 ^ {Jy/z{t))2{Jx{t))2 
for mean values of products of transverse and longitudi- 
nal spins one obtains 



dt 



^jAt) 



f + d(t)FP2(t))Sj,(t) 



Appendix C: Towards purely dissipative steady state 
entanglement using incoherent pump fields 

In the following, we explain how the estimates shown 
in Fig. [3J3 are obtained. We use here a simplified model, 
which has been employed in [Ij to compare the measured 
results to the theoretical predictions. Details regarding 
these fits can be found in the Supplemental Material 
of [JJ . We model the experiment in terms of three atomic 
levels |t>7 = |4,4),, \i)i = |4,3)/ and \h)i = |3,3)/ (|t 
),, = |4,-3),,, \i)ii = \4,-4:)ii and \h)ii = |3,-3),,) 
for the first (second) ensemble. For the timescales consid- 
ered here, the atomic population in other internal states 
is negligible. In order to describe the underlaying physics 
qualitatively, using only a small number of parameters. 



we assume further that F 



and F 



|4,±4)^|/i) 



+ N2{t) (f + d{t)rP2{tf (/i - uf) , viations F 



|/i)^|4,±4) 



■ \h)^\4,±3) 



J- |4,±3)^|/i) — J- out 

where the abbre- 



|4,±4)^|4,±3) 



■ 4,3 



and F 



|4,±3)^|3,±4) 



■ 3,4 
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have been used. 

Atomic transitions are taken into account by introduc- 
ing the collisional rate Fcoi. Since the atomic thermal 
energy is large compared to the level splittings, we as- 
sume the same rate Fcoi for all atomic transitions. Fi- 
nally, we include (t± polarized pump and repump fields, 
which induce resonant transitions with Atti^? ± 1 in the 
first/second ensemble as shown in Fig. [2b. Pump fields 
drive transitions within the manifold of atomic states 
with F = 4 and repump fields transfer atoms from states 
with F = 3 back to F = 4. In this case. 



Tout 
Fin 



rout 
L 



- cob 



repump 



- col- 



where Fl^* is the rate at which atoms leave the two- 
level subsystem due to radiative transitions caused by 
the driving field. Frepump is the rate at which the applied 
repump fields transfer atoms back. Transitions within 
the two-level subsystem occur at the rates 



/i^F and Fp^mp are the driving field and pump field 
induced cooling rates respectively. The heating rate 
caused by the driving field is given by u'^ F. Note 
that the application of pump fields leads to an in- 
creased dephasing rate F. In contrast, repump fields do 
not have an effect on F [3 . We estimate the effect of 
these fields on the dephasing rate by adding 2Fpump to F. 

Fig. (3)3 shows the predicted time evolution of en- 
tanglement for d = 55; 100; 150. The parameters take 
the values used to fit the experimental data measured 
in the absence of additional fields (compare Fig. 2, 
panels a and b in [1 ), F = 0.002ms-\ f = 0.193ms-\ 
Tcoi = 0.002ms~^ and Z = 2.5. The presence of both, 
pump and repump fields is included as described above 



with Fr. 



0.160ms- 



r3,4 
r4,3 



/i' r + Fp 

i^^ F + F, 



■ col 7 
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